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We point out that inflaton decay products acquire plasma masses during the reheating phase fol- 
lowing inflation. The plasma masses may render inflaton decay kinematicaly forbidden, causing 
the temperature to remain frozen for a period at a plateau value. We show that the final reheating 
temperature may be uniquely determined by the inflaton mass, and may not depend on its coupling. 
Our findings have important implications for the thermal production of dangerous relics during re- 
heating (e.g., gravitinos), for extracting bounds on particle physics models of inflation from Cosmic 
Microwave Background anisotropy data, for the production of massive dark matter candidates dur- 
ing reheating, and for models of baryogenesis or leptogensis where massive particles are produced 
during reheating. 

PACS numbers: 98.80.Cq 



^1" 

(N 
l> 
O 
m 
O 

^ ■ 

a: 
i 

<D . 
43 ■ 



X 



I. INTRODUCTION 

At the end of inflation [l| the energy density of the uni- 
verse is locked up in a combination of kinetic energy and 
potential energy of the inflaton field 0, with the bulk of 
the inflaton energy density in the zero-momentum mode 
of the field. Thus, the universe at the end of inflation 
is in a cold, low-entropy state with few degrees of free- 
dom, very much unlike the present hot, high-entropy uni- 
verse. After inflation the frozen inflaton-dominated uni- 
verse must somehow be defrosted and become a high- 
entropy, radiation-dominated universe. 

One path to defrosting the universe after inflation is 
known as "reheating" |2(. The simplest way to envision 
the reheating process is if the comoving energy density in 
the zero mode of the inflaton decays into normal particles 
in a perturbative way. The decay products then scatter 
and thermalize to form a thermal background. 1 

Of particular interest is a quantity known as the reheat 
temperature, denoted as Trh- The reheat temperature 
is properly thought of as the maximum temperature of 
the radiation-dominated universe. It is not necessarily 
the maximum temperature obtained by the universe after 
inflation 0,0, HQ. 

The reheat temperature is defined by assuming an in- 
stantaneous conversion of the energy density in the in- 
flaton field into radiation when the decay width of the 



1 We do not consider here the possible role of nonlinear dynamics 
leading to explosive particle production known as "preheating" 



inflaton energy, T^, is equal to H, the expansion rate of 
the universe. The reheat temperature is calculated quite 
easily Q. After inflation the inflaton field executes co- 
herent oscillations about the minimum of the potential. 
Averaged over several oscillations, the coherent oscilla- 
tion energy density redshifts as matter: p^ oc a -3 , where 
a is the Robertson- Walker scale factor. If we denote as 
pi and ai the total inflaton energy density and the scale 
factor at the onset of coherent oscillations immediately 
after the end of inflation, then the Hubble expansion rate 
as a function of a is (Mpi is the Planck mass) 



H(a) 



'St: pi f ai 



(1) 



Equating H(a) and leads to an expression for ail a. 
Now if we assume that all available coherent energy 
density is instantaneously converted into radiation at 
this value of aj/a, we can define the reheat temper- 
ature by setting the coherent energy density, = 
Pi(ai/a) 3 , equal to the radiation energy density, pn — 
(7r 2 /30)(7*T^ ff , where g* is the effective number of rel- 



ativistic degrees of freedom at temperature Trh. 
result is 
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= 0.2 (™) ' a 1 J 2 ^/M^M P ~i , (2) 

where we have expressed the inflaton decay width as = 

There are various reasons to suspect that the reheating 
temperature is small. For instance, in local supersym- 
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metric theories gravitinos (and other dangerous relics 
like moduli fields) are produced during reheating. Un- 
less reheating is delayed, gravitinos will be overproduced, 
leading to a large undesired entropy production when 
they decay after big-bang nucleosynthesis |£|< The limit 
from gravitino overproduction is T^h ~ 10 9 to 10 10 GeV, 
or even stronger [jj. 

Again, we emphasize that the reheat temperature is 
best regarded as the temperature below which the uni- 
verse expands as a radiation-dominated universe, with 
the scale factor decreasing as g^^^^T^ 1 . In this regard it 
has a limited meaning @>0]. As the scalar field decays 
into light states, the decay products rapidly thermalize 
forming a plasma with temperature T. The latter grows 
until it reaches a maximum value Tmax and then de- 
creases as T oc a~ 3 / 8 down to the temperature Tr#, 
which should not be used as the maximum temperature 
obtained by the universe during reheating. The maxi- 
mum temperature is, in fact, much larger than Trm and 
it is incorrect to assume that the maximum abundance of 
a massive particle species produced after inflation is sup- 
pressed by a factor of exp(— M/Trh)- This has impor- 
tant implications for the idea of superheavy dark matter 
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and baryogenesis 



111. 

The goal of this paper is to present a simple, but rel- 
evant observation that changes the usual picture of the 
temperature evolution during reheating. During the pro- 
cess of reheating the inflaton decay products scatter and 
thermalize to form a thermal background. A thermal- 
ized particle species produced during the first stages of 
reheating acquires a plasma mass m p (T) of the order of 
gT, where g is the typical (gauge) coupling governing 
the particle interactions ^2- This happens because for- 
ward scatterings of fermions do not change the distribu- 
tion functions of particles, but modify their free disper- 
sion relations, producing a plasma mass. The dispersion 
relation can be well-approximated for both scalars and 
fermions by to 2 = k 2 + m 2 (T), where ui and k are the en- 
ergy and the three-momentum of the particle in the ther- 
mal background, respectively. The presence of thermal 
masses imply that the inflaton zero-mode cannot decay 
into light states if its mass is smaller than about gT. 
The decay process is simply kinematically forbidden. 2 

Our observation is that during the reheating stage, 
the inflaton starts decaying and the temperature of the 
plasma rises. If the maximum temperature obtained 
by the universe during reheating, Tmax, is larger than 
about g~ x M^, the inflaton decay channel into light states 
become inaccessible and the decay process stops as soon 
as the temperature has reached a value of the order 
of g M$. Subsequently, expansion cools the plasma, 
lowering the temperature and the corresponding plasma 



masses of the light states. The inflaton is the free to 
decay. However, as soon as this happens, the tempera- 
ture of the plasma rises and the inflaton decay process 
becomes kinematically forbidden again. As a result, one 
expects a prolonged period during which the temperature 
of the plasma is frozen to a plateau value of the order of 
g- x M^. 

Our observation has various implications. First of all, 
let us notice that we do not know the mass of the infla- 
ton field around the minimum of its potential during the 
reheating stage. Indeed, from the recent WMAP cosmic 
microwave background (CMB) anisotropy data EH we 
only have limited informations about that portion of the 
inflaton potential experienced by the inflaton field during 
inflation; we know that it has to be quite flat in order to 
allow a sufficiently long period of exponential growth of 
the scale factor 0,0, EH El- However, we know noth- 
ing about the inflaton mass during reheating, since this 
depends upon a portion of the inflaton potential which 
is not accessible to any observations. This amounts to 
saying that should be regarded as a free parameter. 
Even more, in many inflationary scenarios, e.g., hybrid 
models the reheating dynamics may be determined 
by a scalar field \ different from the inflaton field. (In 
the following, the terminology "mass of the inflaton" will 
be therefore used in a loose way.) 

Suppose that the reheating temperature defined 
in Eq. (J2J), is larger than g~ x M^. This means that when 
the inflaton decay lifetime is of the order of the age of the 
universe, the inflaton field would like to decay, but is not 
allowed to because the plasma masses of the light decay 
products are too large. Only when the energy density 
stored in the inflaton field becomes smaller than about 
P4> ~ (3 _1 -^>) 4 wlu the particles in the plasma have a 
mass smaller than and inflaton can promptly decay. 
Under these circumstances the reheating temperature of 
the universe should be 
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2 This observation was made firs t in in the context of the Affleck- 
Dine baryogenesis scenario ll.'J . 



which is directly related to the inflaton mass and inde- 
pendent of the inflaton decay rate! 

Before concluding the introduction, we note that our 
effect is applicable in situations other than reheating after 
inflation. It would apply, for instance, if the universe is 
ever dominated by a decaying nonrelativistic particle. 

The rest of the paper is organized as follows. In Sec. II 
we analyze in detail the behavior of the temperature dur- 
ing the reheating stage, and in particular we characterize 
the plateau stage both analytically and numerically. Sec- 
tion III is devoted to the study of some applications of our 
findings. We focus on the production of gravitinos dur- 
ing reheating, on the evaluation of the number of e-folds 
after inflation which has recently acquired particular rel- 
evance in order to restrict models of inflation from the 
CMB anisotropy data, and on the production of massive 
particles. Finally, in Sec. IV we present our conclusions. 
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II. REHEATING WITH THERMAL MASSES 

We now discuss the reheating process, assuming that 
the decay products of the inflaton field rapidly thermalize 
and acquire "plasma" masses m p (T) of order gT, where 
g is the coupling constant for a particle in the plasma. 3 

There are two assumptions that deserve elaboration. 
The first aspect is the assumption of "rapid" thermal- 
ization. The timescale for thermalization of the inflaton 
decay products is (ncr) -1 where a is a cross section for 
the scattering of the decay products and n is the number 
density of scatterers. The thermalization is rapid if this 
timescale is short compared to the timescale for energy 
extraction from the inflaton, assumed to be equal to the 
lifetime of the inflaton, T^ 1 . It is reasonable to assume 
that the inflaton is weakly coupled and the inflaton life- 
time is large compared to the thermalization time; hence 
rapid thermalization p| . 

The second important aspect of the assumption is that 
the inflaton decay products have a thermal mass of order 
gT, where g ~ 0.5 is a typical gauge coupling constant. 
One might imagine that the inflaton decays into some 
weakly-interacting particles which then subsequently de- 
cay into "thermal" particles with gauge interactions. But 
in any case, eventually the decay sequence must include 
particles with gauge interactions for which there will be 
a thermal mass. 

To model the effect of plasma masses, let us consider, 
for the moment, a model universe with two components: 
inflaton field energy, p$, and radiation energy density, 
pn, which contains all the light degrees of freedom pro- 
duced after decay. For simplicity, we can think that the 
produced particles in the radiation component all have 
couplings of the same strength. 4 Also, we consider the 
simplest type of decay, that is the decay of the infla- 
ton into scalars. In the case of decay into scalars, the 
only effect of the masses is to modify the phase space of 
the products, while the case of fermions is slightly differ- 
ent, since the scattering amplitude also depends on the 
masses. 

The presence of thermal masses imply that the decay 
width of the inflaton is no longer the zero-temperature 
result r<£ = a^M^, but becomes 



Ta,(T) = a^M^i 
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(4) 

The consequence of this simple fact is that the dynamics 
of reheating drastically changes when the temperature of 
the plasma is such that m p (T) becomes as large as Ma,. 
When m p (T) <C Ma,, the effect is negligible, while the 



3 We expect 3^1, and will assume g = 1/2 for numerical esti- 
mates. 

4 Since the energies are well above the weak scale, even neutrinos 
will interact with gauge-coupling strength. 



decays stop when m p (T) « Ma, since the phase space 
factor goes to zero as T = M^fag ss Ma,. 

With the above assumptions, the Boltzmann equations 
describing the redshift and interchange in the energy den- 
sity among the different components are 



P<P - 
PR 
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- IV(2>* = 



(5) 



where dot denotes time derivative. 

It is clear that the system behaves in such a way that 
T never becomes larger than M^j2g, otherwise the factor 
Ta,(T) would become imaginary. In other words, when T 
reaches this value we have a phase with approximately 
constant T during which the decays are suppressed for 
kinematic reasons. During this phase pa stays constant, 
while p§ decreases like a" 3 . We recall that without 
plasma masses, the behavior of T is very different: im- 
mediately after inflation ends it grows rapidly to Tmax-, 
and then decreases like a~ 3 / 8 until it reaches Trh- At 
this point the <f> field decays completely and the universe 
becomes radiation dominated. 

Taking into account the effect of plasma masses, we 
may have three possibilities: 
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II 
III 



Tmax < M 
Trh < Ma, < Tmax 
Ma, <T rh . 



In case I, the effect of the plasma mass is negligible. In 
case II, after a very short time T grows to Ma,, then stays 
approximately constant for a while, then decreases as 
fl -3/s un tii reheating and the radiation dominated phase 
begins. In case III, again after a very short time, T grows 
to Ma, and after a long phase of constant T, the universe 
directly enters the radiation-dominated phase after the 
time of reheating determined ignoring plasma effects. 

We want now to discriminate, in terms of the funda- 
mental parameters of the inflaton field, the applicable 
case (I, II, or III), and the duration of the constant-T 
phase. First, recall that the maximum temperature ob- 
tained after inflation is given by 6] 

'3\ 2/5 fl5_\ 1/4 1/i (M pi H I " 1 
\ 2tt 3 
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where V is the value of the inflaton potential at the end 
of inflation. The reheating temperature was defined in 
Eq. (J2J). We may now determine the conditions that de- 
termine the operative case in terms of the value of the 
decay constant a. a,: 
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We see that, for realistic values of parameters, it is likely 
that we are in the second or in the third case, i.e., the 
effect is non-negligible. If we put V 1/4 « 10 13 GeV, M 6 w 
10 8 GeV and g* « 10 2 , we obtain 
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II 
III 



a A < 10" 18 
10- 18 < a A <3x 10- ln 
3 x 10- 10 < 
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Next, we may estimate the duration of the constant- 
T phase in cases II and III. We will denote by a/ the 
value of the scale factor at the beginning of the reheating 
phase, and by ajr its value at the end of the constant-T 
phase. 

In case II, of may be estimated by assuming the usual 
scaling of the temperature ignoring plasma mass effects 
during reheating, T oc a -3 / 8 , and finding the value of a 
when T drops below the value Ms/2g. The behavior of 
T is 
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Imposing the condition T/Ms = l/2g to define clf we 
find 



H = (2. 9 ) 8 / 3 ( ~- 
ai x 
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l GeV, 



In terms of number of e-folds, imposing V 1 ^ 4 ss 10 1 
M 6 w 10 9 GeV) we obtain N w 30 + 2/3 ln(c^). 

In case III, the situation is much different than the 
case ignoring plasma effects. In the usual case (without 
plasma masses) the system would enter the radiation- 
dominated era at the time of cf> decay (T^ = H): 
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FIG. 1: The behavior of the temperature during 
without (dashed line) and with (solid line) plasma mass ef- 
fects, for case II: Trh < M^. 
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FIG. 2: The behavior of the temperature during reheating, 
without (dashed line) and with (solid line) plasma mass ef- 
fects, for case III: Trh > Mj,. 
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In our case, though, decays are not possible so long as T 
is larger than M&/2g. So, the <f> energy density continues 
evolving approximately like a~ 3 until p& becomes smaller 
than pr, at which time the <f> can decay without enhanc- 
ing the temperature (and so closing the phase space for 
the decay). So the condition is simply for case III is 



(10) 
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In terms of number of e-folds, imposing again realistic 
values for this case, V 1 ' 4 w 8 x 10 n GeV, M A w 2 x 
10 7 GeV, we obtain N w 14. The two cases reduce to 



the same value in the intermediate case (i.e., the case in 
which Trh — M 6 ). 

Now we want to analyze in detail what happens to 
the system in cases II and III by numerically solving the 
Boltzmann equations. In order to do this it is more con- 
venient to express the Boltzmann equations in terms of 
dimensionless quantities that can absorb the effect of ex- 
pansion of the universe. This may be accomplished with 
the definitions 



R = pro, . 



(12) 



It is also convenient to use the scale factor, rather than 
time, as the independent variable, so we define a variable 
x = aMs- With this choice the system of equations can 
be written as (prime denotes d/dx) 




y/Wx + R. 



$,(13) 
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where the temperature T(x) depends upon R and g», the 
effective number of degrees of freedom in the radiation: 
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It is straightforward to solve the system of equations 
in Eq. l|13|l with initial conditions at x — xi of R(xj) — 
X(xj) = and ^(xj) = $/. It is convenient to express 
P4,{x = xj) in terms of the expansion rate at xi, which 
leads to 



3 M 2 Pl Hj 



The numerical value of is irrelevant. 

We show in Figs. ^ and [21 the solution of the system re- 
spectively in cases II and III. They follow the qualitative 
behavior we described, with the prominent constant-T 
phase. 



III. APPLICATIONS 
A. Thermal production of gravitinos 

The first question we want to address is the produc- 
tion of gravitinos during reheating, taking into account of 
the effect of thermal masses. 5 It is known that the over- 
production of gravitinos represents a major obstacle in 
constructing cosmological models based on supergravity 
0- Gravitinos decay very late and, if they are copiously 
produced during the evolution of the early universe, their 
energetic decay products destroy 4 He and D by photodis- 
sociation, thus jeopardizing the successful nucleosynthe- 
sis predictions @, Oj . As a consequence, the ratio of the 
number density of gravitinos, n 3 / 2 , to the entropy den- 
sity, s, should be smaller than about 



"3/2 



< 10" 



(16) 



for gravitinos with mass of the order of 100 GeV. 

Gravitinos can be produced in the early universe be- 
cause of thermal scatterings in the plasma during the 
stage of reheating after inflation. Usually, to avoid the 
overproduction of gravitinos, one has to require that the 
reheating temperature Trh after inflation is not larger 
than about 10 8 to 10 9 GeV @. In our case, the rele- 
vant parameter is no longer Trh, since the temperature 
is cutoff by the effect of thermal masses. We present here 
an analysis of the thermal generation of gravitinos during 
reheating with a phase of constant temperature. 



5 Here we consider gravitinos, but the results are easily generalized 
to other dangerous light relics. 



Recall the salient aspects of the calculation of the grav- 
itino abundance without thermal masses. The gravitino 
abundance is determined by the Boltzmann equation 
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where (cr^w) oc 1/Mp l is the thermal average of the grav- 
itino annihilation cross section times the M0ller velocity. 
Assuming the actual gravitino density is much less than 
its equilibrium value {n z/2 ) eq = 3g 3/2 C(3)T 3 /47r 2 (g 3/2 
is the number of degrees of freedom of the gravitino), 
the evolution of the comoving gravitino number density 
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(15) (N — a n 3 / 2 ) is quite simple: 
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where c = (35 3/2 C(3)/47r 2 

In the radiation-dominated phase H cx a~ 2 and T cx 
a -1 , so the dominant contribution to -/V 3 / 2 comes from 
small a, corresponding to large T. During reheating H oc 
a" 3 / 2 . If plasma effects are not important T oc a~ 3 / 8 
during reheating, while if plasma effects are important 
T oc const, during reheating. In either case, the dominant 
contribution to iV 3 / 2 comes from large a, corresponding 
to the end of reheating. Therefore we can calculate -/V 3 / 2 
at the end of the reheating era (the beginning of the 
radiation-dominated era), and compare it to the comov- 
ing entropy density N s — a 3 T 3 27r 2 g*/45. The result is 

10-2 ip ( Trh < M * cases r > n ) 

,- 2 M\ l 
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(T RH > M$ caselll). 

(19) 

iaring Eqs. i|16|) and (|19fl . one obtains the bounds 

(10 8 - 10 9 ) GeV > ( T ™ % RH < ¥/ CaSe "ll, II) 
v ; ~ [ Ma, (T rh > Ma, case III). 

(20) 

This calculation illustrates the point that in case III, 
the reheat temperature Trh has no meaning. 



B. Number of e- folds after inflation 

The quality and quantity of observational data has 
reached the point where it is possible to start to place 
meaningful constraints on inflationary models |l4l ITsl 
fill . In the phenomenology of extracting predictions from 
even simple inflation models, one of the significant un- 
certainties is the location of the inflaton corresponding 
to when scales of observational interest crossed the Hub- 
ble radius during inflation. Recent studies of this issue 
[T7L HH have pointed out that a significant factor is the 
uncertainty in the duration of the reheating phase. Lack 
of knowledge of the duration of the reheating results in 
an uncertainty in the number of e-folds of expansion after 
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inflation ends Q{. The uncertainty is usually parameter- 
ized in terms of the reheat temperature, with the uncer- 
tainty in the number of e-folds of inflation depending on 
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As we have stressed, in case III the reheat temperature 
has no meaning; the radiation-dominated era commences 
with T = M^. If case III obtains, then previous formulas 
for the number of e-folds should depend on 



AiV=iln% 
3 V 1 / 4 
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instead of the traditional formula used for AN [2| , AN = 
| IuTrh/V 1 ^ 4 , i.e., T RH ~ y/T^Mpi should be replaced 
by M$. This means that if case III is attained, the num- 
ber of e-folds corresponding to scales of observational in- 
terest is smaller than in the usually adopted case by a 
factor | In ^a^Mpi/M^. 

Proper calculation of the number of e-folds after in- 
flation is crucial in determining the viability of inflation 
models. The change in the number of e-folds in case III 
may be crucial. 



C. Production of massive particles 

Our findings may be relevant for the production of 
massive particles during the reheating stage and, in par- 
ticular, for the production of superheavy dark matter 
(WIMPZILLAS) 0,[I3 and leptogenesis Hj. 

There are many reasons to believe the present mass 
density of the universe is dominated by a weakly inter- 
acting massive particle (WIMP), a fossil relic of the early 
universe. Theoretical ideas and experimental efforts have 
focused mostly on production and detection of thermal 
relics, with mass typically in the range a few GeV to a 
hundred GeV. However, during the transition from the 
end of inflation to the beginning of the radiation phase, 
superheavy and nonthermal particles may be generated. 
If they are stable they may provide a significant contri- 
bution to the total dark matter density of the universe. 

Let us consider a superheavy particle X with mass 
Mx- In this section we will restrict our attention to 
case III for which the final reheating temperature is fixed 
by the inflaton mass, and we consider the case in which 
Mx > M^,. 6 We suppose that the X-particles are pro- 
duced in pairs during the reheating stage by annihilation 
of light states. The corresponding Boltzmann equation 
for the number density nx reads 



dn 



x 



At 



+ 3Hn x = ~{*av) (n x - (n 2 x ) eq ) , (22) 



In the opposite case, the superheavy dark matter may be gener- 
ated directly by the inflaton decay and its abundance will depend 
on the coupling between the inflaton and the X-particles. In such 
a case, the final abundance of the superheavy dark matter will 
be model dependent. 



where (<jav) — ax/M x is the thermal average of the 
annihilation cross section times the M0ller velocity. As- 
suming the actual density nx is much less than its equi- 
librium value (n x )eq = gx{M x T /2Trf/ 2 er Mx ' T (g x is 
the number of degrees of freedom of the A-particles) and 
remembering that dominant contribution to the produc- 
tion comes from end of reheating when the temperature 
is of the order of M^, we can estimate the ratio between 
the number density of X-particles and the entropy den- 
sity at the end of reheating to be 
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corresponding to a present-day abundance of 



!»A^K^^i«, r )(^] e -2M x /M^ (24) 



(23) 



Taking M x (<jav) ~ 1, a moderate hierarchy between the 
inflaton mass and the superheavy dark matter particle 
Mx, Mx/M<p ~ 30, may explain the observed value for 
the dark matter abundance of about 30%. Eq. is 
much different than previous results p|. 

Our findings also have important implications for the 
conjecture that ultra-high energy cosmic rays, above the 
Greisen-Zatsepin-Kuzmin cutoff of the cosmic ray spec- 
trum, may be pro duced in decays of superheavy long- 
living particles |2ll I22I |23| . In order to produce cosmic 
rays of energies larger than about 10 13 GeV, the mass of 
the X-particles must be very large, Mx ~ 10 13 GeV, and 
their lifetime tx cannot be much smaller than the age of 
the Universe, tx ~ 10 10 yr. With the smallest value of 
the lifetime, the observed flux of ultra-high energy cos- 
mic rays will be reproduced with a rather low density of 
X-particles, fix ~ 10~ 12 . The expression Eq. it^Hl sug- 
gests that the X-particles can be produced in the right 
amount by collisions taking place during the reheating 
stage after inflation if the inflaton mass is about a factor 
40 smaller than Mx- 

Let us now discuss the consequences of our results for 
the leptogenesis scenario |20j (even though our findings 
can be easily generalized to any out-of-equilibrium sce- 
nario for the production the baryon asymmetry) where 
the lepton asymmetry (L) is reprocessed into baryon 
number by the anomalous sphaleron transitions |24| . 
Again we will assume case III for which the final reheat- 
ing temperature is fixed by the inflaton mass. 

In the simplest leptogenesis scenario, the lepton asym- 
metry is generated by the out-of-equilibrium decay of a 
massive right-handed Majorana neutrino, whose addition 
to the Standard Model spectrum breaks B — L. 

Let us indicate by tin the number density per comov- 
ing volume of the lightest right-handed neutrino N, the 
one whose final decay (into left-handed leptons and Higgs 
bosons) is responsible for the generation of the lepton 
asymmetry. We can approximate the Boltzmann equa- 



7 



tion for N as 

driN 
dt 



-T N (n N - (n N ] 



(25) 



where IV is the decay rate of N for the processes N — > 
HHl,H£l- Assume again that < Mjv, and that the 
actual density njy is much less than its equilibrium value 



(tin) 



2(M N T/2nfl 2 e- M "' T . 



Since the dominant 



contribution to the production of right-handed neutrinos 
will come from end of reheating when the temperature is 
of the order of Ma,, we can estimate the ratio between the 
number density of TV-particles and the entropy density at 
the end of reheating to be 



tin 
s 




M N /M^ 



(26) 



where in the last expression we have imposed that when 
right-handed neutrinos are produced, their direct decay 
is inefficient, i.e., 



K = 



r 



N 

IT 



T~M<t, 
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1. 



(27) 



The limiting case K ~ 1 would mean that the right- 
handed neutrinos enter into chemical equilibrium as soon 
as they are generated. 

The ratio in Eq. (|26l) remains constant until the right- 
handed neutrinos decay generating a lepton asymmetry 
L = e(n,N/s), where e is the small parameter containing 
the information about the CP-violating phases and the 
loop factors. The corresponding baryon asymmetry is 
B = (28/79)£, assuming only Standard Model degrees 
of freedom, and therefore the final baryon asymmetry is 
bounded to be smaller than 



B = 10" 4 e 



Mn 
Ma 



3/2 



-M N /M$ 



(28) 



For a hierarchical spectrum of right-handed neutrinos, 
it has been shown that that there is a model indepen- 
dent upper bound on the CP asymmetry produced in 
the right-handed neutrino decays, e ^ 3m„ 3 M/v /(87rt; 2 ), 
where m vz is the mass of the heaviest of the left-handed 
neutrinos and v is the Standard Model Higgs vacuum 



expectation value |25| . Therefore, the maximum value 
of the baryon asymmetry in Eq. I128[I is further bounded 
from above by (taking m„ 3 ~ 0.07 eV, the atmospheric 
neutrino mass scale) 



B < 10" 



A I 



N 



10 10 GeV 



Mn_ 

Ma, 



3/2 



-M N /M$ 



(29) 



The requirement that B is larger than 2 x 10 11 implies 
that the ratio Mm I Ma, cannot be larger than about 15. 
IV. CONCLUSIONS 

Reheating after inflation occurs due to particle produc- 
tion by the oscillating inflaton field, and its dynamics is 
very rich. In this paper we have observed that the in- 
flaton decay products acquire plasma masses during the 
reheating phase. The plasma masses may render inflaton 
decay kinematicaly forbidden, causing the temperature 
to remain frozen for a period at a plateau value. This 
happens in any models where the decay rate of the infla- 
ton field r<^ is larger than about M?/Mp;. This condition 
does not seem to be very restrictive. If the condition is 
met, the final reheating temperature is uniquely deter- 
mined by the inflaton mass, and not by its coupling. If 
the reheating dynamics is mainly dominated by a scalar 
field x different from the inflaton, then the final reheating 
temperature may be determined in terms of the mass of 
the x field. An example is if reheating takes place along 
a flat supersymmetric direction whose mass is the soft 
supersymmctry breaking scale fh ~ 10 2 GeV and whose 
couplings to ordinary matter is of order unity. In such a 
case, the effects of plasma blocking are crucial to deter- 
mine the final reheating temperature to be Trh ~ fh. 

We have shown that our results are relevant for the 
thermal production of dangerous relics during reheating, 
for extracting bounds on particle physics models of in- 
flation from Cosmic Microwave Background anisotropy 
data, for the production of massive dark matter candi- 
dates during reheating, and for models of baryogenesis or 
leptogensis where massive particles are produced during 
reheating. 
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